Dynamics of Phase Transitions: The 3D 3-state Potts Model 
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In studies of the QCD deconfmmg phase transition or cross-over by means of heavy ion experi- 
ments, one ought to be concerned about non-equilibrium effects due to heating and cooling of the 
system. In this paper we extend our previous study of Glauber dynamics of 2D Potts models to the 
3D 3-state Potts model, which serves as an effective model for some QCD properties. We investigate 
the linear theory of spinodal decomposition in some detail. It describes the early time evolution 
of the 3D model under a quench from the disordered into the ordered phase well, but fails in 2D. 
Further, the quench leads to competing vacuum domains, which are difficult to equilibrate, even 
in the presence of a small external magnetic field. From our hysteresis study we find, as before, a 
dynamics dominated by spinodal decomposition. There is evidence that some effects survive in the 
case of a cross-over. But the infinite volume extrapolation is difficult to control, even with lattices 
as large as 120 3 . 



> 



o 
in 
o 

o . 

43 ■ 



X 
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I. INTRODUCTION 

In a previous paper p] , which we denote by I in the fol- 
lowing, we investigated the dynamics of 2D Potts models 
in the Glauber universality class. Dynamics in this class, 
model A in the classification of Ref.yK contains local 
Monte Carlo (MC) updating schemes 0] which imitate 
the thermal fluctuations of nature. Here we extend our 
study to the 3D 3-state Potts model with and without 
an external magnetic field. This model gives an effective 
description of the deconfinement transition 0, 0, El 
and by adding an external field |8( , one can represent the 
effect of quark masses. 

A detailed motivation of our study in the context of 
QCD and a compilation of preliminaries such as Potts 
models, spinodal decomposition, Fortuin-Kasteleyn (FK) 
clusters and structure functions, are given in the first two 
sections of I and are not repeated here. Hysteresis cal- 
culations for 2D Potts models are presented in I for a 
number of observables including the (internal) energy, 
properties of FK clusters and structure function. The 
main result is that the dynamics leads for most parame- 
ter values to a spinodal decomposition, which dominates 
the statistical properties of the configurations, indepen- 
dently of whether the equilibrium transition is second 
or first order. Hysteresis cycles and time evolution af- 
ter a cooling quench are presented in I and we find 
also initial evidence (compare Fig. 18 of I) that the sce- 
nario of dynamical spinodal decomposition may survive 
when the equilibrium transition is rendered into a cross- 
over. The latter point is of importance for QCD, because 
lattice simulations and effective models suggest that for 
physical quark masses m u , rrid of the order of 10 MeV 
and m s of the order of 150 MeV neither a chiral nor a 
deconfming transition occurs in the sense that there are 
thermodynamical singularities |l0j| . 



In the linear approximation of spinodal decomposition 
Cahn-Hilliard theory 0, f° r model B predicts an ex- 
ponential growth of the low momentum structure func- 
tion in the initial part of the time evolution. Whether 
such an exponential growth was found or not was used 
by Miller and Ogilvie [3 to determine a critical k c be- 
tween low and high momentum modes, which they re- 
lated to the Debye screening mass. In this paper we 
present the linear approximation, derive its diffusive dif- 
ferential equation, and confront the predicted exponen- 
tial growth with MC data. 

At the present level our results are mainly qualitative. 
Within a model of dynamics one expects universal prop- 
erties, because the equations of motion differ only in their 
constants, but not in their general structure. The proper 
dynamics for Minkowskian relativistic dynamics is ex- 
pected to be hyperbolic. Altogether little is known about 
the kind of effects one may expect to be produced by a 
dynamics. Therefore, we would like to advocate studies 
of different models of dynamics to get a handle on their 
similarities and differences. We begin with the diffusive 
dynamics of model A by the pragmatic reason that it is 
relatively easy to follow the time evolution of correspond- 
ing MC algorithms. The study of hyperbolic dynamics 
within full QCD or just pure non-Abelian gauge theory is 
presently out of reach. However the Polyakov loop model 
of Pisarski and Dumitru 14] allows for simulations in the 
Minkowskian [l5| and it may be interesting to confront 
an investigation of its dynamics with our present results. 

In the next section we fill a gap of the literature by 
presenting the linear theory of spinodal decomposition 
for model A. Our numerical results are compiled in sec- 
tion lllll with quench and hysteresis investigations in sub- 
sections. Summary and conclusions are given in the final 
section ll VI 
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II. SPINODAL DECOMPOSITION IN THE 
LINEAR APPROXIMATION 



Metastable states are non-equilibrium states with life- 
times comparable to the time of observation, while unsta- 
ble states are non-equilibrium evolving states with no ob- 
servable lifetime 0, 0] . When dealing with metastable 
states we consider systems with a first order phase tran- 
sition, while unstable states can be observed for systems 
with first and second order phase transitions. 

Metastable and unstable states are distinguished by 
different types of instability characterizing the early 
stages of a phase separation. Metastable states re- 
quire finite amplitude, localized fluctuations. The sys- 
tem evolves by growing the nuclei or droplets of the new 
phase and the theory of nucleation explains (e.g. Becker- 
Doring theory [lfil fl7l]) the rate of formation of such nu- 
cleating droplets. In contrast, infinitesimal amplitude, 
non-localized fluctuations are sufficient for the initial de- 
cay of unstable states, which is called spinodal decompo- 
sition. 

In mean field theory there is a sharp distinction be- 
tween metastable and unstable states. They are sepa- 
rated by the spinodal curve, which is the locus of points 
inside the coexistence curve for which an appropriate sus- 
ceptibility (for instance the concentration susceptibility 
dc/dfi for binary fluids) diverges. There is no sharp dis- 
tinction in real systems. 

There is a difference in the dynamics for models with 
conserved (model B in the classification of 0) and non- 
conserved (model A) order parameters. Binary alloy or 
binary fluid models are typical examples of models with 
conserved order parameters (the concentration c) . Origi- 
nally IT^ the linear approximation was developed for 
model B. Following we consider here a linear theory 
of spinodal decomposition for model A. Related results 
can also be found in . The general picture is the same 
for model A and model B dynamics. The difference is in 
the quantitative details. 



a statistical treatment of the elementary physical cells. 
The other requirement is that I <C £, where £ is a typi- 
cal correlation length of the system. On the scale I one 
treats the system as being in the equilibrium. The order 
parameter m(r, t) (e.g. the magnetization) has to vary 
smoothly on the length scales larger than I, to satisfy the 
requirement that the system on these length scales is not 
far from equilibrium. 

Since elementary physical cells are equilibrated, their 
equilibrium thermodynamical functions are defined. We 
write down the coarse-grained partition function for the 
system of cells 

/ 

Zi({m}) = Y, c M-E{S}/k B T), (1) 

{S} 

where the sum is over spin configurations {S} and the 
prime indicates the constraint 



(2) 



{S} e cell ft 



By summing in Eq. over all configurations {m} of 
magnetization, one gets the normal partition function. 
The coarse-grained free energy is 



F l {{m}) = -k B ThxZ l ({m}) = J d d rf(m(r)), 



(3) 



where we introduced the non-uniform free energy den- 
sity f(m(r)) and treat the magnetization as a continuous 
function of r. We assume that / is a function of the local 
magnetization and its derivatives and expand it about /o, 
the free energy per spin of a uniform phase. For a cubic 
lattice /(m, Vm, . . .) = fo(m) + k± V 2 m + k 2 (S7m) 2 + . . . 
and 



Fj({m}) ~ / d d r[f a (m) + fciV 2 m + k 2 {Vmf 



(4) 



A. Coarse-grained free energy 

We start by considering the coarse-grained free en- 
ergy It is a functional of the order parameter and 
its minima correspond to the equilibrium states of the 
model. To use the coarse-grained energy for a study of 
non-equilibrium phenomena, one assumes that the model 
is not too far from equilibrium. The requirement is ful- 
filled for systems with slow dynamics on the scale of local 
equilibration times. 

One introduces a spatial mesoscopic length scale I, sep- 
arating the fast dynamics local microscopic regions from 
the slow dynamics macroscopic regions. The system is 
divided into regions (elementary physical cells) of size I, 
which are centered around vector coordinates f defining 
their positions. The length scale I has to be significantly 
larger than the lattice spacing size, I » a, to allow for 



Applying Green's first identity |20j and imposing that 
Vm is zero at the boundary, we obtain 



d d rkiV 2 m 



d d r VfciVm 
d d r ( d ± 



The coarse-grained free energy can be written as 



mm}) 



d d r 



^K(Vm) 2 + / (m) 



1 , dki 
where — K = k 2 : — 

2 dm 



(5) 



(6) 
(7) 



It has the familiar Ginzburg-Landau form. Note that /o 
depends on the choice of the spatial length scale I. 
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B. The linear approximation 

The phenomenological model we use is based on the 
assumption that the rate of displacement of the order 
parameter is linearly proportional to the local thermo- 
dynamic force SF/Sm(f) 



dm 



5F 



dt T Sm{r)' ^ 

The coefficient of the proportionality T, is the response 
coefficient and defines the relaxation time scale of the sys- 
tem. This equation describes a purely dissipative system 
and not one with an internal Hamiltonian dynamics. 

Let us determine the variation of the coarse-grained 
free energy Jfjjl 



6Fi({m}) = / dr 



1 d(Vm) 2 

-K \j- '- 5Vm- 

2 d(Vm) 



dfojm) 
dm 



5m 



Using again Green's first identity and dropping the sur- 
face integral, we write the variation as 

df (m) 



5F t ({m}) 



dr 



-KV 2 m 



The functional derivative is 



SF 



-KV 2 m 



dm 
df (m) 



5m . 



5m(r) dm 
and the equation of motion becomes 



dm 
~dt 



-K V to + — — 

dm 



(10) 



(11) 



(12) 



The Langevin approach 17J gives the same equation, but 
with a noise term added. We are interested in the fluc- 
tuations of the magnetization about some average value 
mo 



(13) 

To solve the differential equation we write it for small 
fluctuations 



m(f,t) = mo + u(r,t). 
ial equat 

du{r, t) 



-KV 2 u(f,t) 



dt 

dfo(m Q ) 
dmn 



(14) 



+ u(f,t) 



d 2 fo(m ) 



dm\ 



and represent u(r,t) as a Fourier series 

u(f, t) = — u(k, t) exp(ifc • r) . 

k 

Using the Fourier series Eq. (|14[) becomes 



(15) 



v 



du(k, t) 
dt 



= -r 



— uj{k) u(k, t) 
dfo(m ) 



exp(ifc • r) 



dm 



where u>(k) = —V ( Kk 



d 2 f Q (m ) 
dm 2 , 



(16) 
(17) 



Multiplying both sides by J^-exp(— ifc' • r) and using the 
Kroneckcr function relation 

V5(k -k')=^ exp (i(k - k') ■ f\ 



we get 



du(k',t) 
dt 



where g(k) = —T 



— uj{k') u{k' , t) = g(k') 



dm a 



The general solution is 



u(k',t) = Ccxp (w(k')t) 



9(k') 

Lu{k') 



(18) 
(19) 

(20) 



Except for the linear term this equation is the same as 
the one produced for model B by Cahn-Hilliard theory 
|llj | and a similar analysis is appropriate. If 



d 2 fo{m ) 
dm 2 ) 



> 



holds, Eq. 1)17(1 implies that the amplitude of any fluctua- 
tion approaches a constant exponentially fast with time. 
But if the second derivative is negative, then one sees 
an exponential growth of the fluctuations for momentum 
modes smaller than the critical value 



1 d 2 f (m ) 



K 



dm 2 



1/2 



C. Equation of motion for the structure factor 

The structure factor (or function) can be measured in 
condensed matter experiments: A system is exposed to 
radiation of a wavelength A and the scattering intensity 
at different angles 6 is recorded. The structure factor 
S(k,t) is proportional to the scattering intensity at the 
angle 8, where k — (4tt/X) sin(i#). In the gauge theory 
studies one may expect that the structure factor reflects 
the production of corresponding momentum gluons. The 
structure function can be written as 



S(k,t) = J {u(r,t)u(r\t)Jexp[ik-(r-r')) drdr' 
= (u(k,t)u(-k,t)) (21) 



To derive the equation of motion for the structure factor 
we take the time derivative of S(k,t) and make use of 
Eq. O 



dS{k,t) 
dT~ 



2u(k,t) u*{k,t)uj(k) 



+ ((u*(k,t)+u(k,t))g(k)) . (22) 
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The average of fluctuations about the average magneti- 
zation has to be zero (u(k, t)) = 0. Thus (|22[) becomes 



dS(k,t) 



Of 



= 2uj(k)S(k,t) 



(23) 



with the solution 



S(k, t) = S(k, t = 0) exp (2u(k)tj . (24) 

We see that it is very similar to the evolution of fluc- 
tuations. Again, if d 2 /o(mo) / drn^ < low momentum 
modes grow exponentially. The value of the critical mo- 
mentum is the same as for the fluctuations. The maxi- 
mum growing mode is the one for which uj(k) takes on 
its maximum, w max = w(fc ma x)- In the linear theory 
femax = according to Eq. (|T7Jl . The values of u(k) de- 
crease with increasing k. Then there exists a k c value, so 
that higher momentum modes k > k c are exponentially 
decaying. The results in the linear approximation are 
the same as for models with conserved order parameter, 
except for the values of Lu(k) and k max . 

A similar structure factor behavior may be obtained 
from the somewhat more accurate Langevin approach 
if we consider only linear terms. This result was first 
obtained for model B by H. Cook according to the re- 
view | 17| . For model A a straightforward algebra yields 



dSjk.t) 
dt 



2u(k)S(k, t) + 2k B TT + 2g(k) . 



(25) 



When the intrinsic dynamics of the system is prevail- 
ing over the extrinsic thermal fluctuations than the right 
hand side is dominated by the first term and the previous 
result for the structure function it^l is recovered. 



D. Validity of the linear approximation 

In the following we assume the existence of two phases, 
each described by a single equilibrium value of an order 
parameter. Potts models belong to a slightly more gen- 
eral group, because there are q = 2,3,... equilibrium 
values in the ordered phase. To trace them in our present 
discussion would lead to a rather cumbersome notation. 
Now the uniform free energy /o at the scale I has two 
minima which correspond to the equilibrium values of 
the order parameter in the two phases: mi and m.2- The 
obvious criterium for the validity of the linear approxima- 
tion is that the mean square value of fluctuations of the 
magnetization are much smaller than the square of the 
magnetization scale, Am — m,2 — mi: {u 2 ) « (Am) 2 . 

Note that the fluctuations u depend on the introduced 
coarse-grained scale. The fluctuations are averaged over 
the cells of volume l d ~ k~ d 

(u 2 {r, t)) = S(0, t) = [ dkS(k, t) ~ k d S(0, t) (26) 



and we can write the condition for the validity of the 
linear approximation as 



k d S(0, t) < (Am) 2 



(27) 



The Langevin approach allows us to estimate the value 
of S(0, t) from the requirement that the first term of the 
right-hand side of i|25|l is the largest (see [ljj for similar 
consideration for model B) and gives 



k B Tk c ! 



k B Tk : 



d-2 



\d 2 f (m )/dm 2 \(Am) 2 K{Amf 



< 1 



(28) 



If we assume that the system is near the critical point 
T c , this equation is not fulfilled. Close to T c the only rel- 
evant length scale is the correlation length £ ~ k^ 1 . The 
analysis of planar interfaces shows that if(Am) 2 /£ 
is essentially a surface tension a in a scaling sense. So 
we can rewrite as 



knTr 



d-l 



<1 



(29) 



But critical point scaling arguments show that the ratio 
(fcsr c )/(cr£ d_1 ) is a universal constant of order unity for 
d < 4. This means that for near-critical quenches the 
linear approximation is not valid and non-linear behavior 
is relevant. 

For quenches far from the transition point T«T C the 
mean-field approach becomes accurate. It allows to eval- 
uate the condition l|28|l of validity of the linear approxi- 
mation. Using mean-field estimates one finds inequality 
to hold for both model A and B 113 : 



I - 



T 



d/2-2 



< 1 



(30) 



where a is the lattice spacing and £o is the range of the 
interaction. For Potts models a/£o = 1- Eq. (|30|l is 
satisfied for d < 4 and T<T C . Note the dependence on 
the number of dimensions in both limits considered. 



III. NUMERICAL RESULTS 

Our Boltzmann weights are exp(— (3E) with the energy 
function 



E : 



(r.r') 



$cr(r,t), a{r' ,t) 



t), (TO 



(31) 



The first sum runs over all nearest neighbor sites r and 
r and a takes the values 1, . . . , q. In this paper we rely 
on symmetric lattices of N = L d spins in d = 2 and d = 3 
dimensions. 

Of the 3D zero magnetic field Potts models the q = 2 
Ising model exhibits a second order phase transition, and 
the transitions are first order for q > 3. The model 
of our interest is the 3D 3-state model, which in zero 
magnetic field has a weak first-order phase transition at 
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C = 0.275283 (6). This transition persists for small val- 
ues of the magnetic field h and has a second order end- 
point at 



C = 0.27469 (1) and h c = 0.000388 (5) 



(32) 



as determined in Ref.|22| (there is a factor of two differ- 
ence between our notation here and this reference). 
The structure factor is (see I) 



S(k,t) 




(33) 



Spinodal decomposition is characterized by an explosive 
growth in the low momentum modes, while the high mo- 
mentum modes relax to their equilibrium values. 

During our simulations the structure functions are av- 
eraged over rotationally equivalent momenta and the no- 
tation i s used to label structure functions of momen- 
tum \k\ = ki where 



k = ^-n 



(34) 



In 2D we recorded the structure function for the modes: 
nf. (1,0) and (0,1), n 2 : (1,1), n 3 : (2,0) and (0,2), n 4 : 
(2,1) and (1,2), n 5 : (2,2). 

In 3D we recorded the modes (including the permuta- 
tions) nr. (1,0,0), n 2 : (1,1,0), n 3 : (1,1,1), ra 4 : (2,0,0), 
n 5 : (2,1,0), ne: (2,1,1), n 7 : (2,2,0), n 8 : (2,2,1) and 
(3,0,0), (3,1,0), mo: (3,1,1), n u : (2,2,2), n 12 : 

(3.2.0) , n 13 : (3,2,1), n u ; (3,2,2), n 15 : (3,3,0), m 6 : 

(3.3.1) , ni 7 : (3,3,2), n 18 : (3,3,3). Note the accidental 
degeneracy in length for n%. 

As in I we measured various properties of FK and ge- 
ometrical clusters. See the definitions given in I. 

For most of our simulations we use the heat bath algo- 
rithm with systematic updating, because we found pre- 
viously that systematic updating is faster than random 
updating. However, for the linear theory of spinodal de- 
composition, one may be interested in the very early time 
development. Then it is advantageous to be able to take 
data within one sweep and we use random updating for 
a few cases. 



A. Quench 

We study the time evolution after a quench from an 
initial to a final temperature 



71 



initial 



71 



final 



(35) 



and are interested in the limit of large lattices, L — > oo. 
To get a statistically meaningful sample, we performed 
between 20 x 20 and 32 x 10 repetitions of each quench. 
Here the first number counts the PCs we used in parallel 
and the second number gives the repetitions performed 
on each PC. 



The large number of repetitions required and the var- 
ious adjustable parameters (Tinitiai, Tfinal and L) make 
a thorough investigation laborious and quite computer 
time consuming. Therefore, we were not able to arrive 
at conclusion for all questions, which deemed interesting 
to us. 



1. Structure Functions 

Miller and Ogilvie investigated the dynamics of 
SU(2) and SU(3) gauge theories after quenching from a 
low to a high physical temperature (corresponding to the 
0mm — > Anax in the spin system). They report a critical 
value k c , so that modes grow (do not grow) exponen- 
tially for k < k c (k > k c ). Using an effective potential 
approach, they were able to relate k c to the Debye screen- 
ing mass. 




100 150 200 250 300 350 400 450 500 
t 



FIG. 1: Structure functions for the 2D Ising model quench 
P = 0.2 -> 0.6 at zero field on an 80 x 80 lattice. 

First we give an example that our 2D Potts model 
data do not support the initial exponential growth pre- 
dicted by the linear approximation. Fig. ^ depicts the 
k > k$ structure function for the 2D Ising model quench 
= 0.2 — > 0.6 at zero field on an 80 x 80 lattice (the 
modes k\ and k 2 are omitted, because their peaks are 
more than four times higher than those of k 3 and k^). 
As these functions do not show an initial exponential in- 
crease, random updating was subsequently used to be 
able to follow the time development within one sweep. 
Figure |21 shows the thus obtained very early time devel- 
opment of the structure functions of S^ for i = 1, 2, 3 
and 4. The increase is linear in t. The curves turn then 
concave, excluding any exponential growth. 

In Fig. [3] we plot our structure functions for a — 
0.2 — > 0.3 quench with sequential updating for the 3D 3- 
state Potts model. Their early time development allows 
for exponential fits of the form 1241) 



Ci(fc) + C 2 (k) exp(2w(fc)i) . 



(36) 
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FIG. 2: Very early time evolution of the 2D Ising model 
quench /3 = 0.2 — > 0.6 at zero field on an 80 x 80 lattice. 
The order in the legend agrees with the order of the curves 
and the lower three curves fall almost on top of one another. 
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FIG. 3: Structure functions for the /3 = 0.2 -> 0.3 3D 3-state 
Potts model quench at zero field on a 40 3 lattice. 
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FIG. 4: The early time evolution for the quench of Fig. [3] 
together with random updating data. The order in the legend 
agrees with the order of the curves. 



indicating an initial, exponential approach towards the 
constant C\(k), reaching for larger times a maximum and 
converting into a decrease on an always concave curve. 
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In Fig. 0] we show the corresponding random and sequen- 
tial updating data (the latter without error bars, not to 
overload the figure). The time between the updating 
schemes is scaled according to i ran d = 1.94 £ soq . 

The fit values for oj(k) are positive for the low values of 
fej, i > 1. Notably, the largest value <x> max is obtained for 
£>3 in contrast to the monotonous decrease of w(fci) which 
the linear theory l|17l) predicts for increasing ki. Already 
our 2D results of Fig. [21 disagreed with the ordering of 
w{ki) predicted by the linear theory, but there it comes to 
no surprise, as the linear theory fails altogether in 2D. In 
3D the &3 mode appears to be the only major exception, 
amazingly already observed as such in pure gauge theory 
[l3j . from a picture which is over-all consistent with the 
predicted decrease of u){k{) for increasing ki. The fits 
become quite unstable when u>(k) approaches zero. For 
even higher ki values negative uj(k) values are obtained, 



FIG. 5: Determination of fc c for the 3D Potts model at zero 
field. 



For sequential updating we have pushed our analysis 
to L 3 lattices as large as L = 110. The results of all our 
uj(k) fits (|36fl are compiled in Fig. [S] versus \n\ 2 defined 
by equation l|34|) . Approximately, we find straight lines 
uj(k) = a + ai \n\ 2 with a negative slope a x and we deter- 
mine the critical momentum k c as the value where u>(k) 
changes its sign. Using l|34f> we find k c ks k§ = 0.351 
(L = 40) and fc c w k w = 0.342 (L = 80). So the finite 
size correction from L = 40 to L = 80 is about 3%. For 
L = 100 and 110 we estimate \n c \ 2 > 31, outside the 
range where we took data. The k c estimate of [H (the 
Polyakov loops substituting for the Potts spins) does not 
consider fnite size corrections, which future studies may 
reveal. 
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2. FK Clusters 

A quench in the temperature changes instantaneously 
the bond probability of FK cluster configurations and ob- 
servables such as the number of clusters, the mean vol- 
ume and surface, the maximum volume, and maximum 
surface change considerably. Then the evolution proceeds 
smoothly. 

The FK Potts clusters are similar to domains in a fer- 
romagnet. For each of the possible magnetizations we 
define the largest cluster. When we quench into the or- 
dered phase one magnetization will eventually take over. 
However, a-priori the system does not know which direc- 
tion this is, because it is prepared in the disordered state 
where all directions are equally probable. Let us address 
the question how the system grows domains. There are 
two alternatives, (i) The system grows a cluster of a par- 
ticular magnetization, while clusters of the other magne- 
tizations decay right from the beginning, (ii) The system 
grows clusters in each direction of magnetization. At a 
later stage smaller clusters coalesce by statistical fluctu- 
ations and one of the directions emerges as dominant. 




FIG. 6: Largest FK clusters of distinct magnetizations for the 
3D 3-state Potts model quench f3 = 0.2 — > 0.3 at zero field on 
various size lattices. 

In Fig. 0we plot the evolution of the largest FK clus- 
ters for the three magnetizations of the 3-state Potts 
model in zero external magnetic field. For different lat- 
tice sizes the plot shows that the system grows clusters of 
each magnetization before one becomes dominant. The 
process of competitions between the largest clusters of 
different magnetization takes longer on the larger lat- 
tices. 

To study the dependence of the speed of the evolution 
on the depth of the quench, we consider different final 
temperatures. In Fig. [7] we compile results for a 40 3 lat- 
tice and find that the growth of the largest clusters of all 
the magnetizations is faster when the system is quenched 
deeper into the ordered phase. 

For gauge theories an equivalent percolation model 




FIG. 7: Largest FK clusters of distinct magnetizations for the 
3D 3-state Potts model quenches from /3 — 0.2 to /3f at zero 
external magnetic field on a 40 3 lattice. 
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FIG. 8: Largest geometrical and FK clusters for the 3D 3- 
state Potts model quenched from /3 = 0.2 to /3f = 0.3 at zero 
field on a 40 3 lattice. 



does not exist and it is problematic to find an appro- 
priate substitute for the FK definition when clusters of 
Polyakov loops are considered. But it is easy to study 
geometrical clusters. Therefore, we compare for the 3D 
3-state Potts model the evolution of geometrical and FK 
clusters in Fig. |SJ The geometrical definition (bonds be- 
tween neighboring spins of the same direction are always 
set) leads to the first of our two scenarios. Geometri- 
cal clusters do not compete. The system starts growing 
one of the domains, while reducing the largest domains 
of the other magnetizations right after the quench. This 
picture is unfavorable for the use of geometrical clusters 
of Polyakov loops in gauge theories. FK clusters are sub- 
sets of geometrical clusters. So, the FK clusters of the fig- 
ure are always smaller than the geometrical. Sufficiently 
far in the ordered phase FK and geometrical domains 
converge in size, while closer to the disordered phase the 
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difference is significant, because geometrical clusters lead 
to artificial groupings of spins into clusters. 
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FIG. 9: The 3-D 3-state Potts model quenched from (3 = 0.2 
to Pf = 0.3 at zero and h — 0.0005 external magnetic fields 
on a 40 3 lattice. 

Next, we apply the external magnetic field h = 0.0005, 
which puts the model slightly above the critical endpoint 
at h c = 0.000388 (05) in the h—T plane, see Eq. (J32J), and 
drive the system through the emerging region of a sharp 
crossover. The external magnetic field sets the direction 
of preferable magnetization. In Fig. [§| we see that the 
growth proceeds in the same manner as for h = with 
FK clusters of all magnetizations growing at the early 
time, only that the preferred domain grows a bit faster 
than before. This indicates that this effect may to some 
extent survive in QCD studies of the crossover region 
and finite size behavior in the L — > oo limit ought to be 
studied, although the physical volume relevant for RHIC 
experiments can only be approximated withing lattice 
gauge theory. 



We build ensembles of at least 160 cycles, grouped in 
bins of 20 cycles each. The error bars are obtained with 
respect to 32 jackknife bins for smaller lattices and 8 for 
larger lattices. The system is driven between /3 m ; n = 0.2 
and /3 max = 0.4. It is equilibrated in the disordered phase 
for at least 80 sweeps at /3 m i n , before being driven through 
a hysteresis cycle. 



1. Latent Heat 

As in I we use the maximum opening of the energy hys- 
teresis as estimator for the latent heat on a finite lattice 
(some energy hysteresis pictures are given in I). Now we 
perform the limit L — > oo for a fixed dynamics defined 
by n'g. The dynamics is at a constant speed in physical 
units, in contrast to the np dynamics used in I, which 
becomes infinitely slow in the limit of large lattices. 



0.14 



g 0.08 




FIG. 10: The 3D 3-state Potts model finite size estimates 
of the latent heat for various n'p dynamics at zero external 
magnetic field. 



B. Hysteresis 

In our 3D hysteresis investigation we look at the first 
order phase transition and rapid crossover region of the 
3-state Potts model. We no longer study the limit of very 
slow dynamics (Eq. (7) of I), but use volume independent 
step-sizes 



A/3' 



2(/3„ 



A. 



(37) 



where Lq is a constant used to set the scale. We chose 
Lq = 20, which gives on a 20 x 20 lattice the same step- 
size as we used in I. The volume independent step-size 
(|37|l makes one sweep a physical time unit for model A 
type fluctuations in nature. It has the nice side effects 
that it allows to simulate larger volume lattices, but ques- 
tions about the survival of effects in the limit of a slow 
dynamics become more difficult to answer. 



Relying on energy hysteresis curves for the 3D 3-state 
Potts model at zero external magnetic field, we plot in 
Fig. Un| the lattice size dependence of their maximum 
openings. A plateau is reached for L > 30 for the n' g = 1 
dynamics and for L > 40 for the large values of n'p = 4 
(slower dynamics). The error bars given in Fig. 1101 arc 
purely statistical and do not take into account a possible 
systematic error, which comes from the L independent 
constant step-size of the n' g dynamics |25j. 

As plateaus for Ae7(n^,L) are reached at L > 40, 
it is save to take the L = 80 latent heat values as 
the our final answers. They are listed in the follow- 
ing, n'p = 1: 0.1249(7), n' a = 2: 0.1128(7), n' a = 4: 
0.1017(5), n' 3 = 8: 0.0932 (6), and n' fi = 8: 0.0895(6). 
They are much lar ger than the equilibrium estimate 
Ae7 = 0.05354(17) 23]. The slower speed of the tran- 
sition at higher n'p lowers the values, but in the limit 
n'p -K»we still expect a dynamical latent heat which is 
larger than the equilibrium result. Namely, the L —> oo 
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limit of the dynamics of I provides a lower bound and 
a dynamically generated latent heat was observed there 
for 2D Potts models, including the Ising model. In 3D 
one may want to modify the dynamics which gives a 
lower bound from A/3 ~ l/V of I to A/3 ~ 1/L 2 , be- 
cause the slowing down at the critical point of the Ising 
model scales with L z and z w 2 (see Ref. |24| for es- 
timates of z). Fitting our present Aej(n^) data to the 
form Ael(n' p ) = AeJ+a (n' /3 ) b gives Ae7 = 0.057 (9) with 
an exponent b = —0.298 (55). The rather large error bars 
prevent a conclusion about whether there is a dynami- 
cally generated latent heat in 3D or not. A change of 
the dynamics and/or a push towards larger values is 
required. 




0.2 0.220.240.260.28 0.3 0.320.340.360.38 0.4 
P 



2. Structure Functions 

The structure function behavior is similar to the 2D 
case with pronounced peaks in the cooling half-cycle. In 
Fig. ^2 we plot Sk 1 for the n'g = 1 dynamics of the 3D 
3-state Potts model in zero magnetic field. In Fig. 1121 we 
show Sfej in the external magnetic field h — 0.0005, which 
is slightly higher than the critical end-point field. Recall 
that ki depends on the lattice size L through Eq. (|34|l . 
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FIG. 11: Hysteresis of the structure function on vari- 

ous L 3 lattices for the 3D 3-state Potts model in zero external 
magnetic field. 

Arrows indicate the flow with increasing /?. At (3 — 0.4 
the hysteresis turns around and the lower curves describe 
the /3 max — * Anin approach. For the smaller lattices (most 
clearly for L — 20) they exhibit small peaks around the 
transition temperature, which are similar in size to peaks 
one finds for equilibrated configurations. For the larger 
lattices these peaks are not seen. Their returning curves 
come in too high, because their hysteresis are at (3 = 0.4 
still far out of equilibrium. 

As expected, the /3 m i n — > /3 m ax peaks are smaller in the 
presence of the external magnetic field. While for h = 
they increase over the entire range of our lattices, we find 



FIG. 12: Hysteresis of the structure function Sk ± (ft) on vari- 
ous L 3 lattices for the 3D 3-state Potts model in the external 
magnetic field h = 0.0005. 



for h = 0.0005 an initial increase, which turns around be- 
tween L = 40 and 60. However, the subsequent decrease 
appears to be slow compared to non-critical behavior for 
which the structure functions falls off ~ l/V in our 
normalization. Possibly, the behavior is similar to the 
equilibrium behavior at a second order phase transition 
point, where the fall-off is ~ l/L x and < x < 3. 
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FIG. 13: Structure function Ski maxima versus lattice size 
for the n'fj — 1 dynamics of the 3D 3-state Potts without 
external magnetic field (h = 0) and in the cross-over region 
(h = 0.0005). For L — 20 and 30 the data points are slightly 
shifted, so that it becomes visible that there are results for 
both h = and h = 0.0005. 

Including lattices as large as 120 3 , we plot in Fig. ITTA 
our estimates of the structure function maxima for Sk ± 
(h = and h = 0.0005) together with fits for the larger 
lattices. Unfortunately, the precision of the data does not 
allow to determine the exponent x in a fit of the form 



C'liiax 



= di + a 2 L x 



(38) 
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accurately. Choosing the exponent by hand, x = +1 for 
h = and x = — 1 for /i = 0.0005, we find a satisfac- 
tory goodness-of-fit Q (for the definition of Q see, e.g., 
Ref. |23) hi each case. Using the non-critical x = — 3 
(instead of a; = —1) to fit the h = 0.0005 data, gives 
the too small goodness-of-fit Q = 0.0062. Of course, this 
cannot exclude that the behavior turns to non-critical for 
even larger lattices. In this connection it is puzzling that 
the value x = +1 used to fit the h = data cannot be 
asymptotic. In our normalization the Sk structure func- 
tions are bounded from above by Sk < 1, implying x < 
for any true asymptotic behavior. Apparently lattices of 
size L > 120 are needed. With checkerboard updating 
they could be handled on supercomputers . 
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FIG. 14: The structure function Sk± hysteresis for the 3D 3- 
state Potts model in the h = 0.0005 external magnetic field on 
a 40 3 lattice for simulations with the dynamics indicated 
by the extension to the lattice size. The inlay enlarges the 
equilibrium peak together with the slowest dynamics np = 16 
data (P is mapped on 0.21 +2(/3- 0.26) and S kl on 4S kl + 
0.0001). 

For a 40 3 lattice we plot in Fig.^Jthe structure func- 
tion Ski for the 3D 3-state Potts model in the h = 0.0005 
field for different up up to 16 together with equilibrium 
data. The approach to equilibrium with the decrease 
of the speed of the dynamics (increase of n'p) appears 
to be faster than in 2D. Possibly a greater connected- 
ness of domains allows in higher dimensional models for 
a faster evolution of macroscopic structures. The heat- 
ing (/3 max — ► 0min) peak develops and reaches the equi- 
librium peak, as we reduce the speed. The cooling peak 
decreases, but it would need much slower simulations to 
reach the equilibrium values. Thus within a range of 
speeds the cooling peak stays strong and spinodal de- 
composition is the dominant scenario. For h = the 
figure looks very similar, but for an increase of all peak 
values by a factor of about two. The explanation of the 
difference between the heating and cooling branches pre- 
sented for 2D models in I remains valid. Driven from 
the disordered to the ordered phase, the system freezes 
in domains of different magnetization, which are slow to 



evolve. For a fast dynamics it does not have enough time 
to equilibrate in the ordered phase. The decrease of the 
speed allows then the system to equilibrate in the ordered 
phase and the heating peak emerges. 
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FIG. 15: Hysteresis for the structure functions for the 3D 
3-state Potts model on an 80 3 lattice in zero external field. 
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FIG. 16: Hysteresis for the structure functions for the 3D 3- 
state Potts model on an 80 3 lattice in the h — 0.0005 external 
field. 



In Fig. E| and ^| we plot for an L — 80 lattice the 
lowest momentum (k±, feg) structure functions for zero 
and non-zero external magnetic field. In the absence of 
the external magnetic field we observe that the peaks 
are very pronounced for k\ to /C3, and much smaller for 
k > ki. In the presence of the external magnetic field all 
peaks are smaller and the difference between lower and 
higher modes is less distinct. We still observe a qual- 
itative difference between the k < k$ and the k > k± 
modes. Another observation is that the peaks are less 
broad for h = 0.0005 than for h = 0, because the system 
equilibrates for h > easier in the ordered phase. 
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3. FK Clusters 

In I we identified the maximum connected surface of 
FK clusters as an interesting quantity. In the transition 
region it exhibits a pronounced peak related to percola- 
tion. The 3D analysis |l8j| of this quantity has remained 
limited to rather small lattices (L < 30), because of a 
rapid slowing down with lattice size of the algorithm 
identifying the maximum connected surface. In contrast 
to 2D deviations from equilibrium appear relatively weak. 
A reason may be that the extra dimension adds new de- 
grees of freedom to the surface. For the crossover region 
(h = 0.0005) the approach to equilibrium is illustrated 
in Fig. ll7l bv comparing the dynamics at different speeds 
(including a 'super-fast' n'p — 1/4 choice) with equili- 
brated data. Notably the figure is practically identical 
for the h = transition, reiterating that percolation hap- 
pens for a crossover in quite the same way as for a proper 
transition |27| . 
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FIG. 17: Hysteresis for the largest connected cluster surface 
fimax for the 3D 3-state Potts model in the external magnetic 
field h = 0.0005 on a 20 3 lattice. The hysteresis flow is indi- 
cated by the arrow. The n'p values are given the extensions 
to the lattice size and e stands for equilibrium. 



IV. SUMMARY AND CONCLUSIONS 



ever, the ordering of the growth coefficients is not always 
in agreement with Eq. I|17ll . Further, the linear approx- 
imation fails for near critical quenches (as was known 
before [Ti|l and, altogether, in 2D. 

Following the structure function evolution over ex- 
tended times, their pronounced peaks are the most no- 
ticeable feature. It is certainly a challenge to base a the- 
ory of spinodal decomposition on their finite size scaling 
instead of using the linear approximation. Unfortunately, 
very large lattices appear to be needed to extrapolate the 
infinite volume limit for the structure function maxima. 
Within our hysteresis investigation we pushed their anal- 
ysis to systems as large as 120 3 . But, for zero external 
magnetic field we find an increase, which cannot yet be 
asymptotic, because it violates an upper bound. Ap- 
plying a small external magnetic field (h = 0.0005) and 
quenching through the cross-over region results in an in- 
crease with L for small lattices (L < 40) and a slower 
than non-critical fall off on our larger [L > 50) lattices 
(Fig. 113(1 . If the last trend continues towards L — > oo, 
signals of spinodal decomposition would survive the con- 
tinuum limit of an analogue QCD scenario. 

A quench into the ordered phase leads to a competition 
of FK cluster domains. For the proper transition (h = 0) 
this leads to a divergence of the equilibration time in the 
limit L — > oo, an effect well-known in condensed mat- 
ter physics |23| ■ For a cross-over the L — > oo scenario is 
less scary, but our Fig. [§] shows that one may still expect 
a substantial time delay before reaching equilibrium. In 
the RHIC experiment one performs a near-critical quench 
into the quark-gluon plasma phase (assuming the phase 
is reached). For a near-critical quench FK clusters are 
initially small. In our Fig. [7| the largest one covers only 
about 10% the system. So, a number of vacuum domains 
are expected and the effects from the dynamics are man- 
ifested on their surfaces. Inside the domains the system 
may equilibrate fast, but the system as a whole cannot 
reach equilibrium easily. Unfortunately, a concise def- 
inition of FK clusters is not available for lattice gauge 
theory. Therefore, one may turn towards investigating 
directly the SU(3) gluonic energy density for influences 
of the supposedly underlying vacuum structure. First 
results will be published in Ref. p9j . 



In our simulations of the model A dynamics of the 
3D 3-state Potts model on L 3 lattices, we find spinodal 
decomposition to be the dominant scenario of phase con- 
version. Under a quench from the disordered into the 
ordered phase, we observe an early time development of 
structure functions Sfc(i), which is in over-all agreement 
with the exponential growth predicted by the linear the- 
ory of spinodal decomposition for k < k c (Fig. [SJ . How- 
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